Ground state properties of a trapped few-Boson system under rotation 
— beyond the "lowest Landau level" approximation — 
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We consider a harmonically trapped few-Boson system under rotation and investigate the ground 
state properties beyond the usual "lowest Landau level" approximation by using exact diagonaliza- 
tions in a restricted Hilbert subspace. We find that both the effective interaction energy and density 
distribution are strongly affected by the two-body interaction strength. 
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Following the experimental realization of Bose-Einstein 
condensation of alkali- metal atoms Q, there has been 
much attention attached to the behavior of these sys- 
tems under rotation both experimentally and theoret- 
ically. Mattews et al. Q created vortices in a two- 
component system, and Madison et al. || studied the 
rotations in a stirred one-component Bose-Einstein con- 
densate. Theoretical studies have mainly focused on the 
Thomas-Fermi regime of strong interactions or on 
the limit of weak interactions between the atoms [p|4l5[ . 
In the latter case, most research works have applied the 
"lowest Landau Level (LLL)" approximation, in which 
all the bosons are in single-particle orbitals with radial 
quantum number n — and angular momentum m hav- 
ing the same sign as the total angular momentum. The 
intermediate regime of moderate interaction, however, is 
rarely concerned [ |l6| , |l7| |. 

In this Brief Report, we would like to study ground 
state properties of a trapped rotating Bose-Einstein con- 
densate with an arbitrary two-body interaction strength 
beyond the LLL approximation. To attack the problem, 
we calculate low-lying states of a few-Boson system by us- 
ing the exact diagonalization method. We assume that 
our numerical computation is qualitatively applicable for 
a rotating Bose-Einstein condensate though the number 
of bosons is limited to a small value of N = 16. We also 
assume that the cloud of atoms rotates about some axis, 
and that the system is in its ground state with respect 
to this axis, which implies that our problem is essentially 
two dimensional. 

The trapped Bose-Einstein condensate, comprised of 
N repulsively interacting alkali atoms with mass M 
and s-wave scattering length a SCl obeys the many-body 
Hamiltonian 
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where the energy and length are given throughout in 
scaled harmonic oscillator units Hlu and a = y/h/Mui, 
respectively. The two-body contact interaction strength 
is characterized by a dimensionless coupling constant 



g = a sc /a. As we shall see, the problem has three differ- 
ent interesting regimes parametrized by the product gN: 
the Thomas- Fermi regime gN 3> 1, the moderate inter- 
action regime gN ~ 1, and the weak interaction regime 
gN « 1. To diagonalize the Hamiltonian, Eq. (1), the 
single-particle states of the Fock state are chosen to be 
eigenstates of the two-dimensional harmonic oscillator, 
i.e., 



= 2ra- 
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where Ljl n '(r 2 ) is the associated Laguerre polynomial, n 
and m are the radial and angular momentum quantum 
number, respectively. Previous studies in weak interac- 
tion limit usually take single particle states with n = 
and m — 0, +1, +2, +3,..., which corresponds to the LLL 
approximation. Here we set up Fock states by extending 
to n ^ 2 and —2 ^ m ^ 4 |lq] , and sample over the full 
Hilbert space with a fixed number of bosons. From this 
sampling, only those Fock states with a given total orbital 
angular momentum L and a configuration energy (corre- 
sponding to the sum of occupied single particle energies) 
less than or equal to a specified cutoff energy E c are in- 
cluded (see figure 1 below for example). The purpose was 
to select only the most important Fock states from the 
full basis, thereby reducing the matrix dimension to an 
acceptable size d < 10 5 . Once the active Fock states are 
constructed, we calculate the matrix elements and sub- 
sequently diag onalize the matrix by using the Davidson 
algorithm jljj , which is very efficient to solve the eigen- 
value problem with a large and sparse matrix. 

It is important to point out that the above diagonal- 
ization scheme only in principle yields an exact solution 
of the many-body problem. For reasons of numerical fea- 
sibility it is necessary to truncate the set of basis func- 
tions to be used in the diagonalization. One then has to 
make sure that convergence of the ground state energy 
is reached with respect to the cutoff. As the required 
matrix size increases rapidly with gN and L, computa- 
tional expenses severely restrict the calculations to only 
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the smallest systems at not too large values of g. Thus, 
with increasing boson number or g, the results become 
less accurate due to the restricted number of basis states 
that can be included in the calculations. Figure 1 shows 
the convergence of the ground state energy as a function 
of the cutoff energy E c for a system of N = 10 bosons at 
L = 10. The lowest possible Fock state has all ten bosons 
distributed in the lowest Landau level, and the corre- 
sponding configuration energy equals to N+L = 20. This 
means that for the ground state energy with different cut- 
off energies displayed in figure 1 all excitations up to an 
energy E c — 20 are included. It is readily seen that both 
the case with a moderate interaction strength gN = 1 
(figure la) and a relatively strong interaction strength 
gN = 3 (figure lb) show convergence at E c > 30. The 
smaller g is, the more rapidly the convergene reaches. 
Therefore, we conclude that our numerical results shown 
below is very accurate for a weak or moderate interac- 
tion strength, and also qualitatively accurate in the rela- 
tively strong interaction regime. In what follows, we take 
E c = N + L + 10. 

We now analyze the influence of interaction strength g 
on ground state energies as a function of the total angular 
momentum, which may be written as 

E gs (L) = (N + L) + gV int ,L, (4) 

where Vmt.L is introduced as a scaled effective interaction 
energy. Figure 2 shows the L dependence of Vint.L of a 
system of TV = 16 bosons for g = 0.006, 0.1 and 0.3. For 
comparison, the result with LLL approximation is also 
displayed. Note that the case of g = 0.006 is closer to ex- 
perimental values for 87 Rb p0| ], i.e., a sc ~ 100a = 5.29 
nm and a — 1.25 x 10" 4 cm. For all three interaction 
strengths, the interaction energy Vi n t,L simply decreases 
linearly or near linearly with increasing the angular mo- 
mentum. At N = L, a kink appears in slope. This is a 
hint of condensation into a vortex state: in macroscopic 
superfluids, the state for L — N would have a condensate 
of unit angular momentum and would be lower in energy 
than neighboring ground states. In the weak interaction 
limit (gN = 0.096 <C 1), our results only show slightly 
deviation from that with LLL approximation. However, 
with increasing the interaction strength, the interaction 
energy Vi n t,L drops significantly for each angular momen- 
tum. This implies a spatial reconfiguration of bosons 
since they should push away from each other to reduce 
the strong interaction as g increases. 

We next investigate the density distribution of bosons 
for various interaction strengths as shown in figure 3. 
In general, the density is found to become decreasingly 
small and broad as g increases. For all angular momenta, 
the shape of curves for the weak and relatively strong in- 
teraction strength differs largely: not only the maximum 
density of the former case is nearly two times larger than 
that of the latter one, but the maximum position is also 
shifted strongly. Moreover, when the angular momentum 
is closer to L = N, a maximum surrounded by an out 



ring of lower density appears at the r = for a relatively 
strong interaction strength g = 0.3 (see figure 3c). 

To further examine the structure of the ground state, 
we calculate the conditional probability distributions 
(CPDs), which measure the density correlation among 
bosons Ipi]] . Unlike the usual density distribution, which 
is circularly symmetric under the rotation invariant con- 
finement, the CPD reflects an intrinsic density distribu- 
tion of bosons Q . And thus it is indeed an observable 
quantity in experiments instead of the circularly sym- 
metric density distribution. This is best illustrated for 
a condensed boson gas, where the circularly symmetry 
is broken spontaneously, and CPD is then related to the 
measurement of boson density p2| . Note also that the 
CPD has been extensively used in the study of electron 
correlations in doubly excited helium-like atoms |25| ] and 
analysis of the formation of Wigner molecule in quantum 
dots p3). We define the CPD for finding one boson at r 
given that another is at ro as 

<* GS | E ^r-r^ro-r,) | * GS ) 
V(r ' = (iV-l)WslE^o-r,)|vI/ GS ) ' (5) 

3 

where I^gs) represents the ground state. In figure 4, we 
show a plot of CPDs for some selected angular momenta 
and for a system of N = 16 bosons. With increasing the 
two-body interaction strength, the bosons spread more 
widely in x-y plane as expected, in line with density dis- 
tributions shown in figure 3. On the other hand, the evo- 
lution of producing a single vortex state at N = L does 
not change significantly even as g is varied over several 
(~ 2) orders of magnitude. 

In conclusion, we have studied the ground state of a 
repulsively interacting Bose-Einstein condensate with a 
nonvanishing angular momentum. We found that all 
physical quantities, especially the effective interaction en- 
ergy, are strongly affected by a dimensionless interaction 
strength g. We realize that our result is based on a few- 
Boson system, and its validity to a Bose-Einstein conden- 
sate should be further checked by more rigorous analytic 
and numerical treatments, such as the quantum Monte- 
Carlo simulation. We also note that although the so-far 
published experiments have mainly studied systems that 
do not satisfy the condition gN ~ 1 or gN <C 1, systems 
that do satisfy this condition are accessible with current 
experimental techniques |25[| . 
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Figures Captions 

Fig.l. Ground state energy of a single vortex state 
versus the cutoff energy for a system of N — 10 bosons 
at the moderate (a) and relatively strong (b) interaction 
strengths. The energy is measured in units of Hlj. Both 
the two cases show convergence at E c > N + L + 10. 

Fig. 2. Vi n t^L (in units of Huj) as a function of angular 
momentum for different two-body interaction strength. 
In the LLL approximation, it is described precisely by 
an algebraic expression Vi n t,L — N(N — L/2 — 1) for 
2 < L «C N ||Jl|]. 

Fig. 3. Density distribution of N = 16 bosons for var- 
ious angular momenta. The unit of length is the oscil- 
lator length a = \J h/Muj. The short dashed, dotted, 
dash dotted lines correspond to g = 0.006, 0.1, and 0.3, 
respectively. The results with LLL approximation is also 
delineated by a solid line for comparison. 

Fig. 4. Selected conditional probability distributions 
for a system of N = 16 bosons. The value of x, y in each 
subplot ranges from —3.0a to +3.0a. ro = (0, 1.0a), and 
other choice of r does not alter our results qualitatively. 
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